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I. 



INTRODUCTION 



In an earlier work Bohannan [2] developed a mathemati- 
cal model for a system of stored items which are periodically 
reworked to improve their reliability. That model yields 
an expression, for the reliability of an item following 
rework, in the form of a series. More recently, Bishop [1] 
applied the discrete control theory to systems encountered 
in economics and operations research rather than the more 
typical applications in electrical systems. 

In this thesis the mathematical model for the system of 
stored items which are periodically reworked is developed 
in terms of control theory, and the solution obtained for 
item reliability has a closed form that is more suitable for 
further studies in this area. Important characteristics 
such as system stability and the steady state value of the 
reliability are directly derived from the control model 
rather than obtained by inspection of simulated values. 

The reliability of the stored items will depend upon 
the reliability at the time of acquisition together with 
the storage environment and length of time the items are 
stored. When periodic rework of items is done, other varia- 
bles that can affect the reliability of the item are the 
effectiveness of the rework, the rate of rework, the rate 
of acquisition of new items and the rate and policy of 
expenditure of items for use, obsolescence, or even training 
purposes . 
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Such a system might include a stock of ordnance which 
is acquired, stored, and periodically reworked, but not 
expended except for war time use, or alternatively expended 
in training with replacement by new items . 

On the other hand, items such as big missiles are not 
likely to be expended, and thus expenditure and replacement 
cannot be considered as a control variable for the relia- 
bility. Thus, unless due to other factors, the expenditure 
of such items would be avoided in favor of rework of the 
items . 

In this thesis we will study only the case where there 
is no expenditure and replacement, leaving reliability over 
time to be maintained by rework. We will also consider 
the case where we have several levels of rework and want 
to find the optimal sequence of levels over successive 
reworks in order to optimize a given performance measure. 

A general rework model is developed in Chapter II 
which relates the reliability of an item following rework 
to its reliability following the last rework. An equiva- 
lent model is also developed relating the reliabilities 
of an item immediately before consecutive reworks. These 
models may be solved for any rework cycle given the initial 
reliability at acquisition. 

In Chapter III the concepts of control theory are 
reviewed and applications to control of the reliability of 
items in inventory are suggested and structured. 
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When we have random errors in our measure of the item 



reliability or when the rework process introduces random- 
ness in the item reliability, then we have to make our 
judgments based on estimated values. In Chapter IV, a 
Kalman Filter is used with our control model to accomplish 
this estimation. An examination of the final value and the 
stability of our system is followed by an extension to 
more than one dimension for the case where we have systems 
composed of several subsystems in series, parallel or a 
combination of both. 

In Chapter V we shall extend our work to the case where 
we can decide among several levels of rework. Using dynamic 
programming, we will show how to derive the optimal rework 
policy for this case in accordance with several performance 
measures, namely, to obtain a desired reliability, to mini- 
mize the costs of several reworks, to minimize the time 
to achieve a desired reliability, or to satisfy a combina- 
tion of all these performance criteria. Conclusions and 
recommendations for further study are offered in Chapter VI. 



i 

i 
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II. DEVELOPMENT OF THE REWORK MODEL 



In this chapter the general rework model is developed 
which will relate either the reliability of an item follow- 
ing a rework or the reliability of an item immediately 
before a rework to the same reliability for the last rework 
and the effectiveness of the rework process. A solution 
is found for these two reliabilities at any instant given 
the initial or acquisition reliability of the item. 

A. THE EFFECTIVENESS OF REWORK 

One form of rework mechanism would raise an item's 
reliability to a certain level which is independent of the 
item's reliability prior to rework, as might be the case 
when components or parts are replaced. Another way, 
developed by Bohannan [2] and adopted in this thesis, is 
for the rework mechanism to achieve an increase in relia- 
bility which is proportional to both the item unreliability 
before rework and the effectiveness of the rework mechanism. 
This type of rework mechanism might exist where major assem- 
blies or subassemblies are tested and repaired rather than 
replaced. Denoting the item reliability immediately before 
the rework by R, the reliability following the rework by 
R s , we can define the effectiveness of the rework process 
a as 

R - R 

d = , 0 < a £ 1 (1) 
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The interpretation of the rework effectiveness a is that 
the greater its value, the more effective the rework pro- 
cess or in other words the item unreliability, 1-R, will 
be reduced by an amount proportional to a since from 
Equation (1) we can write 



Under these circumstances, it may be possible for such a 
rework mechanism to improve the item reliability to where 
it is "better than new", "as good as new" or "not as good 
as new". Because of deterioration of reliability during 
storage, we note that in this last case, the reliability 
might continue to decrease until the items need to be 
replaced, rather than reworked. 

B. THE REWORK MODEL 

From the definition of effectiveness of the rework 
process , given by expression (1) , the reliability of 
an item following the rework, R , is 



1-R 



(1-R) - a ( 1-R) . 



s 



R 



s 



R + a ( 1-R) 



or 



R 



s. 



(1-a) R + a . 



( 2 ) 
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Let R(t) represent the reliability of an item of age t 
when there is no rework. We assume that all failures are 
random, suggesting (i) that there are no early failures or 
that some form of "burn in" has been used to eliminate 
early age failures, and (ii) that because the items are 
in storage and not in an operating environment, either 
there are no wearout failures or the time to occurrence of 
wearout is much longer than the projected time until the 
next scheduled rework. Thus assuming that the reliability 
R(t) is an exponential function, 



R (t) = e _(a+bt) 



(3) 



for t >_ 0 . Since the initial reliability R q , for t = 0, 
is 



R q = R ( 0 ) = e 



we have 



R ( t) = R e bt (4) 

o 

for t _> 0 and b > 0. The value of the parameter b depends 
on the nature of the item stored and on the storage environ- 
ment. 

If the system contains N items, which are reworked at 
a constant rate p, then the period T to "turn over" the 
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inventory is 



We begin by looking at new items that are going to be 
reworked for the first time. The duration of the rework 
process is T/N and the "age" of a new item selected for 

rework under a First In, First Out (FIFO) policy will thus 

T . . • 

be T - jy, and its reliability R immediately prior to its 

first rework will be 

R (T-|) = R c e S (5) 

If the item is going to be periodically reworked, we 

can rewrite the above equation for the reliability just 

prior to the (k+l)st rework. This will occur at time 
T 

t = (k+l)T- jj, and we have 

T 

R(kT + T - 7j) = R s (kT) e 



When k = 0 in this equation, we get 



R( T -|) 



R s (0) e 



-b(T-^) 



and from Equation (5) we have then 



R 



o 



B s<°> ' 



12 



or in other words the initial or acquisition reliability 
may be considered as the "reliability following the 0 tb 
rework", or as the initial value for the reliability 
following the rework, R s (k). Evidently there is no rework 
at the time of acquisition and this equivalence is used 
only for purpose of coherence. The reliability just prior 
to the first rework is then only defined for values of k 
equal to or greater than one. 

The process is presented schematically in Figure 1. 

Assuming that N or the rework rate are relatively large, 

T. - 5 ~ T and 
N 



R (kT+T) = R s (kT) e 

For simplicity we will denote R(kT) by R(k) , and the 
reliability immediately before the (k+1) st rework, R(k+1), 
can then be rewritten as a function of the reliability 
following the kth rework R s (k) : 

R (k+1) = R s (k) e" bT (6) 

This expresses the reliability of an item immediately 
before the (k+l)th rework recursively, as an exponential 
function with initial value R g (kT) which is the reliability 
of the item at the beginning of the kth period. 

Rewriting expression (2) for the reliability following 
a rework with this notation, we have 
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Turn over period T 

— 




R(k) R s (k) 



i T/N ; 

f* * 

1 1 3* tine 

R(k+1) Rg (k+D 



Rework Tine 



FIGURE 1. Reliabilities during the kth period 
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( 1-a) R(k) + a . 



(7) 



R s (k) = 

Combining the above equations, we have the expression 
for the reliability after the (k+l)th rework in terms of 
the reliability after the kth rework, 

—hi 1 

R (k+1) = (1-a) R (k) e + a, k = 0,1,2, ... 

(3) 

We may also write the reliability immediately before 
the (k+l)th rework in terms of the reliability just prior 
to the kth rework, 

R(k+1) = R(k) e~ bT + a (l-R(k) ) e" bT , k=l,2,... (9) 

When using this equation we have to consider as initial 

condition the reliability immediately before the first 

-bT 

rework, R(l) = R e 

o 

C. THE SOLUTION TO THE REWORK MODEL 

Equations (8) and (9) are first order difference equa- 
tions, and we can solve them by using the z transform tech- 
nique, described in Appendix A, with the advantage that we 
will obtain a closed form solution instead of a solution 
in the form of summation. 

In Appendix A we show that taking the z transform of 
both sides of equation (8), we get 



L 
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zR g ( z ) - zR q = (1-a) R s (z) e _bT + j , (10) 

where R (z) means the z transform of R (k) . 

b b 

Rearranged, this gives 

(z - (1-a) e bT )R s (z) = zR Q + p 



or 



V z) * 



z - (1-a) e 



-bT R o + 



az 



(z - ( 1-a) bT ) (z-1) 



which when expanded in partial fractions yields 



R s (z) 



R + 



a 



,, > -bT o z-1 , . , . -bT r ,, . -bT, 

z-(l-a)e l-(l-a)e [z-(l-a)e ] 



a 



r , , , -bT , , 

[ ( 1-a ) e -1] 



The solution of this equation is obtained by taking 
its inverse z transformed, as developed in Appendix A, 
yielding 



o n \ . k -bkT 

R_(k) = (1-a) e R - 

o O 



a 



-bT 

(l-a)e -1 



rn , n .k -bkT, 
[l-(l-a) e ] 



;id 



Using the same method for Equation (9) , we take the 
z transform to get 



— hT — hT 7 

zR ( z ) - zR (1) = ( 1-a) e R(z)+ae , (12) 



16 



or 



[z-(l-a)e~ bT ]R(z) = zR ( 1) + . 

Solving for R(z), we get 

-bT 

R(Z) = =E5 R(1) + — ^ =bT“ 

z-e (1-a) (z-1) [z-(l-a)e Di ] 



which expanded in partial fractions yields 



R(z) = 



z-e bT (l-a) 



•R(l) + 



ae 



-bT 



(z 1} [ 1- ( 1-a) e bT ] 



-bT,,, , -bT, 

ae [ ( 1-a) e ] 



-bT -bT 

[z- (1-a) e Di ] [ ( 1-a ) e -i; 



Taking the inverse z transform, we find the solution in the 
time domain to be 



n m \ / -i , k -bkT„ , , ae 

R(k) = (1-a) e R(l) - 



-bT 



, , , -bT , 

(l-a)e -1 



,, ,, ,k -bkT. 

Ll-(l-a) e 



(13) 



Equation (13) gives the reliability of an item immediately 
before the rework at time t = kT for a given initial con- 
dition R(l), the reliability of the item just before the 
first rework. Equation (11) gives the reliability of the 
item after the rework at time t = kT for a given initial 
condition R q , the reliability of the item at the time of 
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acquisition. Comparing these two equations we notice that 



been expected from the Expression (6) that relates these 
two reliabilities. 

As an example, values of R (k) and R(k) are shown in 
Table (I) for various values of R , a, b and T. As can be 



seen in Table I, for T = 400 and R q = 0.8 we have, for 

a = 0.8, the case "not as good as new"; for a = 0.569 the 

case "as good as new"; and for a = 0.7, an improvement of 

item reliability constituting the "better than new" case. 

The value of a that leads to the case "as good as new" can 

be found by setting the reliability following the first 

rework R (1) equal to the initial reliability R . Thus 
s o 

from equation (11) we have 



-bT 

they are different only by a factor of e , as might have 



o 




(l-a)e bT R(0) 




R ( 0 ) , 



(l-a)e -1 



which, rearranged, gives 



a 



(l-a)e bT ] R ( 0 ) 



or 



a - cxR ( 0 ) e 



-bT 



-bT 

R ( 0 ) (1-e Di ) . 



Solving for a yields 
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k 

1 

2 

3 

4 

k 

1 

2 

3 

4 

k 

1 

2 

3 

4 



TABLE I 



Examples of reliability immediately before rework 
R(k) and reliability following rework R s (k). 
Initial reliability = R q 
R ework Effectiveness = a 
Deterioration parameter = b 
Turnover period = T 



R =0.8 T = 400 

o 



a = 


0.2 


a = 


0.569 


a = 


0.7 


R(k) 


R s (k) 


R(k) 


R s (k) 


R(k) 


R s (k) 


0.536 


0.629 


0.536 


0.800 


0.536 


0.861 


0.422 


0.537 


0.536 


0.800 


0.577 


0.873 


0.360 


0.488 


0.536 


0.800 


0.585 


0.876 


0.327 


0.462 


0.536 


0.800 


0.587 


0.876 




R = 
o 


0.8 


a = 0.7 






T 


= 50 


T = 


100 


T = 


250 


R(k) 


R s (k) 


R(k) 


R s (k) 


R(k) 


R s (k) 


0.761 


0.928 


0.724 


0.917 


0.623 


0.887 


0.883 


0.965 


0.830 


0.949 


0.691 


0.907 


0.918 


0.975 


0.859 


0.958 


0.701 


0.912 


0.928 


0.978 


0.866 


0.960 


0.710 


0.913 




II 


0.95 


T = 100 






a 


= 0.07 


a : 


= 0.8 


a = 


0.9 


R(k) 


R s (k) 


R(k) 


R s (k) 


R(k) 


R s ( k) 


0.860 


0.958 


0.860 


0.972 


0.860 


0.986 


0.867 


0.960 


0.979 


0.976 


0.892 


0.989 


0.869 


0.961 


0.883 


0.977 


0.895 


0.990 


0.869 


0.961 


0.884 


0.977 


0.895 


0.990 



b = 0.001 for all cases 
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a 



(14) 



R (0) 

1 - 



(1-e 
R (0 ) e 



-bT 

-bT 



) 



D. THE CONTROL VARIABLE 

As pointed out by Bohannan [2], with these values a 
steady state in reliability is reached by the fourth rework 
and the influence of the initial reliability is lost. Thus 
in order to rework an item to achieve a desired reliability, 
we have to choose values for the rework effectiveness a, 
or the turnover period T or both, so that the steady state 
is not below the desired value. The turnover period can 
be shortened by increasing the rework rate p or by reducing 
the inventory size N. Another way to have a higher relia- 
bility would be by improving the storage environment, which 
would reduce the reliability deterioration rate b, if this 
is possible. However, the inventory size N and the deteriora- 
tion rate b are not generally considered part of a rework 
policy. In our formulation only the rework rate p and the 
rework effectiveness a could be changed within a rework 
policy. 

In this thesis we choose to work with a fixed rework 
rate p and to use the rework effectiveness a to control the 
item reliability. This appears to be the usual practice in 
preventive maintenance involving a fixed schedule and 
several levels of maintenance or rework (corresponding to 
several values of a) . The use of the rework rate p to 
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control the system makes the system non linear, involving 
equations that are time variant. This can be solved by 
another area of the control theory which while not covered 
in this thesis, can be done in a further work. 

In this chapter we derived a rework model for the case 
when we have a constant rework rate. We solved our dis- 
crete recursive equations for the reliability before the 
rework, R(k) , and after the rework, R (k) , by the z trans- 
form technique. The result is in a closed form rather than 
in the form of a summation, making this solution more suit- 
able for further analysis as will be done in Chapter IV. 
Before this analysis, we will first derive in the next 
chapter this same closed form solution by applying the 
control theory approach. 
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III. DEVELOPMENT OF THE CONTROL MODEL 



In Chapter II we developed the rework model and found 
the system of first order difference equations relating 
the reliability immediately before the rework R(k) to the 
reliability after the rework R (k) : 

b 



R (k+l) 



R s (k) 



e 



-bT 

r 



(15) 



and 



R s (k) = a (l-R(k) ) + R(k) • 



(16) 



We then solved this system using the z transform technique, 
obtaining an expression for the reliability after any 
rework as a function of the initial reliability, R , 



,, v , . > k -bkT^ 

Rg(k) = (1-a) e R Q - 



a 



, . . -bT . 

( 1-a) e -1 



n , n ,k -bkT, 
[1- (1-a) e ] , 



(17) 



and an expression for the reliability immediately before any 
rework as a function of this reliability before the first 
rework : 



R(k) 



(1-a) 



k -bkT 
e 



R(l) 



ae 



-bT 



(1-a) 



-bT 

e 



-1 



[1- ( 1-a) 



k -bkT. 
e ] 

(18) 
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In this chapter we will introduce a control theory 
approach to the rework problem, and then derive these same 
expressions in terms of control theory. This will provide 
a basis for later solving the problem of obtaining an 
optimal sequence of rework levels. 

A. THE CONTROL THEORY APPROACH 

In this section we will review, in a simplified form, 
the main aspects of control theory that are of interest 
to our rework model. 

The first fundamental concept underlying a control 
system is that of a dynamic model, that is a model des- 
cribing the relationships among the relevant variables and 
parameters of our system. In the model, the variables are 
allowed to change with time in a deterministic way, and 
the system behavior depends not only on their values at any 
instant, but also on their past values and the rate of 
change of these variables with time. The equations derived 
in Chapter II constitutes such a dynamic model of our system. 

The second fundamental concept we will need is that of 
a closed-loop control system , i.e. , one in which the output 
has a direct effect upon the control action, as shown 
schematically in Figure 2. The output and the reference 
input are compared in the controller and we say that the 
output is fed back to the controller. The difference between 
the reference signal and the feedback signal, called the 
error signal or input signal to the plant, is then used to 
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FIGURE 2. Block Diagram of a Control System 
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change the behavior of our plant or process with the goal 
of bringing the output of the system to a value that gives 
a zero error. 

The "plant" is any physical or abstract object to be 
controlled. Sometimes we feedback a measure of our output: 
for this case, a box representing the measurement device 
would be drawn in the feedback branch of the block diagram 
of Figure 2. In summary, the term "closed loop" implies 
the use of feedback action in order to reduce system error. 

An open loop control system is a control system in which 
the output has no effect upon the control action, that is, 
the output is neither measured nor fed back for comparison 
with the reference. An advantage of the closed-loop control 
system is that the use of feedback makes the system output 
relatively insensitive to external disturbances and internal 
variations in system parameters, since a correction action 
is taken anytime the error is not zero. 

An useful tool in control theory is a signal flow graph . 
This is a pictorial representation of a set of simultaneous 
algebraic equations in which each variable is represented 
by a graphical symbol called a node, and the dependencies 
between pairs of variables are represented by directed 
branches drawn between pairs of nodes. These dependencies 
between two variables are called transfer functions, or 
gains and are defined as the ratio of the incoming variable 
to the variable at the end of the branch. Transfer functions 
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are used to label the branches , and indicate that a multi- 
plication operation is done upon the value entering the 
branch in the arrow direction, delivering a new branch value 
to the node where the branches terminates. The nodes are 
also summing devices which sum all values arriving by the 
way of incoming branches. As an example of this, Figure 3 
shows the signal flow graph for a hypothetical system 
defined by the set of equations 



x^ = bx^ + dx^ + fu 



x 



2 



ax 



1 



x 3 = ix 2 + ex 4 + gu 



x 



4 



cx 



3 



Y = hx 4 



Given the signal flow graph of a system we often wish 
to write the overall transfer function from an input or 
reference value to the output. This can be obtained 
directly from the signal flow graph by inspection or by 
use of Mason's Theorem. This theorem, which we will use 
later, expresses the transfer function in terms of various 
loop gains and the parallel gains from input to output, and 
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g 




FIGURE 3. Example of a signal flow graph diagram 
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states that the transfer function from input u to a response 
Y is 



T 



Y 

U 



l hh 



where the terms are defined as: 

(1) A is the determinant of the feedback configuration 
and is calculated from the equation 



A = 1 - k. + L.L. - L L L + ..., 

£ 3 L k l L mno 

where the L. (or L, , etc.) are loop gains all the way 
J K 

around a feedback loop in the system. Thus means the 

sum of all loop gains. The next term J' L^L, is the sum 

of all products of pairs of different loop gains - e.g., 

L 1 L 3 anC ^ S ° The P r ime the summation means we use 

only the products for pairs of gains of non-touching loops. 

In other words, is included only if loop 1 does not 

touch loop 3 (two loops touch if they have at least one 

node in common) . Likewise, j ' L L L is the sum of all 

'-■mno 

products three at a time, where again each of the three 
loops does not touch the other two . 

( 2 ) is a path gain of the ith forward path from 
input U to output Y (a path which contains no loops) . 

(3) A^ is the system determinant A after we have 
excluded all loops which touch the path. 
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The Mason's Theorem is better understood by an example. 
Suppose our system is described by the signal flow graph 
of Figure 3 and we want to determine the overall transfer 
function Y/U. The loop gains are 



= ab , L 2 = ce , L 3 = aicd . 



Hence 



[L. = ab + ce + aicd . 

L 3 

The loop L 3 touches both and but and are 
non- touching. Hence 

F'L k L 2 = aicd , 

and there are no three non-touching loops. We can then 
write the determinant A 

A = 1 - JLj + , 



or 



A = 1 - (ab + ce + aicd) + aicd . 



The direct paths from U to Y are 
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f a i c h 



P 



1 



and 



P 2 = g c h . 



Since P^ touches all loops, 



Path P 2 does not touch , hence 



A 2 = 1 - = 1 - ab . 



Substituting into Mason's formula, we have 



T = - 

U 



P 1 A 1 + P 2 A 2 



or 



T _ faich + g c h ( 1 - ab) 

1 - (ab + ce + aicd) + aicd 

Summarizing, we can say that Mason's theorem gives a simple, 
fast procedure for writing a desired transfer function 
directly from the signal flow diagram. The response of 
our system can then be immediately derived from this trans- 
fer function as proportional to the input , 
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Y = T U , 

and this is called the solution to our control model, 
i.e., an expression for the response of our system as a 
function of any input to the system and a known (derived by 
Mason's Rule) overall transfer function. 

B. THE RELIABILITY CONTROL MODEL 

Our concept of a system is not limited to physical 
systems. The concept can be applied to abstract, dynamic 
phenomena such as those encountered in economics and opera- 
tions research. Also, feedback control systems are not 
limited to the field of engineering but can be of particu- 
lar interest to the manager, public official, operations 
researcher, biologist, and design engineer. In this section 
we will structure the reliability model in control theory 
terms. The plant of the control theory approach will be 
the inventory, where the (decreasing) reliability of items 
is our variable of interest, or state variable . The con- 
troller will be substituted by a decision-making function 
that will decide upon a level of maintenance which is optimal 
in the sense of minimizing a certain performance measure. 

We will defer optimization until Chapter VI, and for the 
present will assume a single level of rework, as defined 
in Section D of Chapter II, corresponding to a certain 
effectiveness of the rework a. 
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Our dynamic model is described by expressions (15) and 



( 16 ) 



R (k+1) 



R s (k) e 



-bT 



k = 0,1/2 



and 



V k) 



=>(1 - R (k) ) + R(k) , k = 1,2,3, ..., 



and the signal flow graph can be drawn in several ways, one 
of which is shown in Figure 4. At the second node from 
the left the present item reliability, whose measured value 
is fed back to this point, is compared with the reference 
value, R d , that is made here equal to 1.0 due to our 
definition of rework effectiveness (this will be better 
clarified later). The difference, when not zero, produces 
an automatic decision to rework the item, and after the 
rework, this difference becomes multiplied by a since from 
expression (16) we have 



This difference will be denoted by U (k) , and called 
the input to our plant. Physically the item is returned 
to the inventory at this point and the reliability of the 
item is given by the sum 



R s (k) - R (k) 



a(l - R (k) ) 



U (k) 



(19) 
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FIGURE 4. Reliability Flow Graph for the Rework Model 
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U (k) + R (k) . 



R s (k) 



This sum is performed at the fourth node. While in the 



inventory the item reliability will deteriorate by a 



factor e 



-bT 



given by expression (15) and at the end of a 



period of length T we have the item reliability R(k+1) 
immediately after the (k+l)th rework represented at the 
fifth node by R(k+1). The relation between R(k) and R(k+1) 
is a function called a unit delay of time, represented by 
a box in the flow graph, with the meaning that the output 
of our system R(k+1) is one period T ahead of the value 
fed back, R(k) , used to compute it. 

The variable Y(k) is used to emphasize the fact that 
R(k+1) is the output of our control system. Branches with 
a gain equal to 1.0 are used when there are variables that are 
just renamed or going to be operated on in another mode. 

Taking the z transform on both sides of Equations (15) 
and (16) by applying the properties given in Appendix A, 
we have 



R(z) 




( 20 ) 



and 




a R ( z ) ) + R ( z ) , 

z-1 



( 21 ) 



and the signal flow graph can be drawn as in Figure 5 . 
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Reliability Signal Flow Graph to the Control 
Model in the z Domain 
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The response of the system, Y(z), can be now obtained 
by applying Mason's theorem as shown earlier. Since we 
can consider that we have here two input values, namely 
R(l) and R Q , this yields the overall transfer function 
with respect to the initial condition R(l) as 



Y(z) 

R(l) 



-bT -l -bT .. . -bT 

(1 - ae z +e z-(l-a)e 



and with respect to the reference value R^, 



-bT. z . -1 

YU) ( FT> z 

r d ' l- ae - b V 1 + e - bT 



-bT 

ae z 



(z-l) (z-(l-a)e bT ) 



The response of the system due to the reference value 
that here is made equal to 1.0 is then 



Y (z) 



_ -bT 
R~ae z 



(z-l) (z- (1-a) e _bT ) 



and the response of the system due to the initial condition 
is given by 



Y (z) 



RUJ z 

-bT 

(z- (1-a) e Di ) 



Using the superposition principle for differential 
equations, the total response is 



Y (z) 



-bT 
ae z 



R(l) z 



-bT -bT 

(z-l) (z- (1-a) e ) (z-(l-a)e ) 



( 22 ) 
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This is the same as the z transform of Equation (13) 
in Chapter II, and the inverse z transform of the total 
response then gives the same Equation (13) for the relia- 
bility of an item immediately before a rework. This is 



Y (k) = R (k) = (l-a) k e _bkT R(l) — 



-bT 



-bT 

(l-a)e -1 



ri ,k -bkT, 

[1- (1-a) e ] , 



k 1,2,3, ... , 



where the argument (k+1) was changed to k. In this case 
the z transform had to be applied because the unit delay 
present in the system, does not provide a linear function 
necessary for the application of Mason's theorem. In the 
z domain, however, the unit delay becomes a very simple 
function and the solution may be found easily. 

If the reference input R^ is not equal to 1.0 but is 
the desired value of the item reliability that we want 
R s (k) to match after some reworks then the definition of 
rework effectiveness would become 

R s (k) - R (k) 

“ = R d - R (k) * 

Here, if the reliability immediately before a rework was 
greater than the value of the desired reliability, the 

I 

value of a would be negative. To avoid this problem we 
might use a decision rule that could be expressed as a 
function M(R, R^) defined as 
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0 if 



R d - R (k) < 0 



M 



< 



( 23 ) 



a 



if Rp - R (k) > 0 



with the meaning that no rework has to be done when the 
item reliability R(k) is greater than the desired value 
Rp. Our model would then become 



This approach, however, would introduce a non-linearity 
into our system that could not be eliminated by the z 
transform procedure, and the solution of our problem would 
probably have to be obtained by computer. Hopefully we 
can handle this problem and find an optimal sequence of 
reworks that makes the value of R g (k) to match the desired 
value R d by using the optimization process that will be 
developed in Chapter V. 

Summarizing, we have structured the item reliability 
model in terms of control theory for periodic reworks with 
effectiveness of rework a as defined by equation (1) . We 
are now able to apply the tools of control theory to find: 
(i) the limits of stability of our system as a function 
of the rework period T, (ii) the stable final value of the 
reliability as the number of reworks increase, and the use 
of a Kalman filter. This will be shown in the next chapter. 



R s (k) = R (k) + M(R d - R (k) ) 



( 24 ) 
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IV. APPLICATIONS OF THE CONTROL THEORY MODEL 



In this chapter we will show how the results of our 
control model can be applied to study the behavior of the 
item reliability. First, we will study the stability of 
our system, and then the final value reliability when t 
increases. After introducing randomness to the system, 
we will show the use of a Kalman Filter to estimate the 
value of the reliability that is fed back. 

A. SYSTEM STABILITY 

A linear control system is stable if the output even- 
tually comes back to an equilibrium state when the system 
is subjected to a disturbance. In our reliability system 
we are interested in knowing if there is any range of values 
of the rework effectiveness factor a for which the system 
becomes unstable, i.e., the reliability might continue to 
decrease until the items need to be replaced rather than 
reworked. 

From control theory we know that the characteristic 
equation of a control system is equal to the denominator of 
the overall transfer function plus 1.0 and that for the 
system to be stable the roots of the characteristic 
equation, in the z domain, have to be inside the unity circle 
[10]. For our reliability system we have from Equation (21) 
the characteristic equation in the z domain as 
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(z-1) [z- (l-a)e" bT ] +1 = 0. 



Rearranging our characteristic equation, we have 



z 2 - [1+ (1-a) e -bT ] z + [ 1+ ( 1-a) e~ bT ] = 0 , 



with complex roots at 



-bT * 

= ± j j J [1+ (1-a) e -bT ] 2 -4 [1+ (1-a) e -bT ] 



The absolute value of these complex roots is given by 



Here , 



or 



/l, n I 2 -2bT " . 
= yj[(l-a) e -1] 



z < 1 



z <1 



when 



1 r .2 -2bT . . . . 

2 [d-a) e -1] < 1 



or 
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This gives a second order inequality in a 



2 



2a + 1 - 3e 



2bT 



< 0 



a 



that when solved shows us that for a stable system, the 
value of a must be within the range 



that agrees well with the value of a in Table I for which 
the item is "as good as new" after each rework. 

Outside these limits or when a is less than 0.56, (since 
a by definition is less than unity) , the system is unstable. 
Here, the item will have its reliability deteriorating with 
time until the item needs to be replaced rather than 
reworked, i.e., for a < 0.56 eventually R -+■ 0 . 

B. THE FINAL VALUE IN STABLE SYSTEMS 

When the system is stable its response tends to a 
constant value as the time increase whenever there is no 




(25) 



For the example of Chapter II with b = 0.001 and 



T = 100, we get 



0.56 < a < 2.56 
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disturbance input. To find this steady state value we 
will apply the final value theorem, derived in Appendix A, 
that gives us the limit of the system response when time 
becomes very large. The theorem is stated as 



lim Y (k) = lim [(z-l)Y(z)] . 

k-»- « z->-l 



For our overall transfer function, this results in 



lim R (k) 

k+co 



az 



lim [— 
z^l z-(l-a)e 



, zR ( 0 ) (z— 1) , 

-bT .. , -bT J 

z- (1-a) e 



a 



1- ( 1-a) e 



-bT 



(26) 



Substituting selected values from Table 1, we get 



a = 0 . 7 


T = 100 


lim 

t“M» 


R s (k) 


= 0.961 


a = 0.8 


T = 100 


lim 

t+oo 


V k) 


= 0.977 


r- 

• 

o 

II 

3 


T = 400 


lim 
t ^°° 


R s (k) 


= 0.876 



and these limits agree well with the values taken from 
the table for the fourth rework. 

As we can clearly conclude from expression (26) the 
steady state value of the reliability after the rework R S ( OT ) 
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depends only on the rework effectiveness a and on the turn- 
over period T for a given degradation parameter b, and is 
thus independent of the initial value of the reliability 
as was intuitively observed in Chapter II. 

C. ERROR COMPENSATION USING A KALMAN FILTER 

In this section we will show how a Kalman Filter can 
be applied to our model in order to compensate for errors 
in the measurement of the item reliability. Again, measure- 
ment here does not imply a physical action. We suppose the 
nature of the stored item is such that it receives frequent 
diagnostic checks which yield an estimate of its reliability, 
and such "measurements" are subject to random error. 

In the model developed in Chapter III we have shown how 
the item reliability is fed back so that a decision can be 
made about doing a rework so that the item's reliability 
can approach a desired value. In practice we have to 
measure the item reliability or estimate it by some way in 
order to feed back this observed value. This measurement 
process, however, contributes to the variance within the 
system. Other sources of variance in the process include 
variability in the item environment and in subsystems 
reliabilities due to repairing or substitution of parts. 

We call these random errors as noisy or random input to 
our plant, since we can model these errors as random 
variables that are input somewhere in the plant as will 
see later. 
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If we do not want our control system following these 
random inputs excessively it would be a good policy to 
smooth the reliability values over time and then feed back 
a predicted value upon which we made our decisions. 

A useful/ smoothing device is a Kalman Filter, and a 
block diagram of a general system with a Kalman filter is 
shown in Figure 6. Other filters could be used here but 
with a discrete system we chose the Kalman filter: it is 

the optimum recursive filter in the sense that it minimizes 
the variance of the estimator error. The random input at 
the plant in Figure 6 is designated by the random variable 
w(k) and represents the error introduced by the environment 
and rework process. The random error in the measurement 
is designated by the random variable v(k) . 

Making the assumption that these noises are additive 
expected values, we can restate our system equations (15) 
and (16) as 



R(k+1) = e bT R s (k) + w (k) 



and 



R s (k) = ( 1-a) R (k) + aR D , 

and substituting the second equation into the first, we 
have 
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w(k) 




FIGURE 6. Block diagram of a general control 
system with Kalman Filter 
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R (k+1) 



[ (l-a)R(k)+aR D ] + w(k) . 



= e 



-bT 



Rearranging and grouping terms we have 



R (k+1) = e" bT R(k) + e -bT U(k) + w(k) (27) 



where 



U (k) = a[R D - R(k) ] 

as before and R^ = 1. 

If we introduce the Kalman Filter in our system we 
are feeding back an estiamted value R(k) instead of the 
actual value R(k) and the input to the plant becomes 

U (k) = [R D - R (k) ] (23) 

The measurement process can be modelled by the equation 

Y (k) = H R (k) + V (k) (29) 

where H can be any function describing the measurement, and 
v(k) is the random variable in the measurement process. 

We need the following additional assumptions to make 
the application of a Kalman Filter valid: 

(a) The measurement error v(k) and the random process 

input w(k) has zero mean and are uncorrelated between 
states, with variances S and Q respectively, or 



46 



E[v(k) v ( j ) ] 



S ^ > j — 0,1,2, 



E[w(k)w(j) ] 



Q ^ > 3 — 0,1,2, 




(b) The initial state is a random variable which has 
known mean E[R(0)] = R q and variance M. 

(c) The estimator is characterized by the linear 
relationship 



R (k | k) = R(kjk-l) + G(k)[Y(k) - HR(k|k-l)] , k = 0,1,2,. 



R(k|k-1) is the optimal prediction of R(k) given 
observations at times up to and including k-1, and 
G(k) is the gain of the Kalman Filter at each step. 

(d) The random input and initial state are uncorrelated. 

(e) The random errors in the plant process or random 
input w(k) and the random error in the measurement 
process v(k) are uncorrelated or 



where: 



R(k|k) is the optimal estimate of R(k) given 
observations at times up to and including k, 



E [w (k) v ( j ) ] 



0 



j ,k = 0,1,2, . 
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The functioning of the Kalman filter is by predicting 
the value of the item reliability R(k[k-1) using the esti- 
mated value at the end of the (k-l)th period, multiplying 
this predicted value by the function H, so that this new 
value can be directly compared with the measured value. 

The filter gets the correction term [ Y (k) -HR (k | k-1) ] 
that can be conveniently weighted by the gain G(k) to 
correct the predicted value, and obtains a new estimated 
value for the kth period. 

Notice that if we did not use a predicted value, 
what is equivalent to say that R(k|k-1) = R(k-ljk-l) = R(k-l), 
we would have exponential smoothing and the gain G(k) 
would then be the smoothing constant. With the predicted 
value any trend is recursively incorporated in the esti- 
mated value. 

A 

The best prediction R(k|k-1) that can be made is 
clearly by using our model (Equation (27)), 

A ^ , -hT 

R (k | k-l) = e R (k-1 1 k-1) + e U(k-1) , (30) 

and from assumption (c) above, the equation for the esti- 

/\ 

mated value R(k|k) is 



A A A A 

R (k) = R (k | k) = R (k [ k-1) + G (k) [Y (k) -HR (k | k-1) ] . 

(31) 

With these two equations and the assumptions above it 
is shown [9] that the sequence of gains that minimizes the 
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variance of the estimator error P(k|k) is given by the 
recursive equations 



G (k) = P (k | k-1) H [H 2 P (k | k-1) +S ] _1 , (32) 

P (k | k) = [1 - HG ( k ) ] P (k | k-1) , (33) 



and 



P (k | k-1) = (e~ bT ) 2 P (k-1 | k-1) + Q. (34) 

These are initialized by the value 

P ( 0 |-1) = M . (35) 

From these equations we can see that the Kalman Filter 
gain G(k) does not depend on the measurement values Y (k) , 
and can thus be computed in advance and stored for later use 
when processing the measurements as they become available. 

We can see also that the gain at time k, G(k), is inversely 
realted to the variance of measurement error S — the more 
uncertainty in the measurements, as reflected by a larger 
S, the smaller G(k) will be, and the less the correction 

A 

term [Y (k) -HR (k | k-1) ] in Equation (31) will be weighted in 
determining the next estimate. The random input error Q 
also affects the gain G(k), and as the uncertainty in the 
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prediction of P(k|k-1) increases and the more uncertainty 
we have in our model, the larger will be G(k) . For large 

A 

values of M or no confidence in the initial guess R(0|-1) 

/v 

the filter makes R(0|0) equals to the first measurement 

Y(0) since G(0) is equal to 1.0. When the initial guess is 

better the filter makes a weighted combination of the guess 
/\ 

R(0|-1) and the first measurement Y(0), as one can conclude 
from Equation (31) . 

A flow graph in the time domain of our system is shown 
in Figure 7. The lower part of the flow graph is the 
Kalman Filter, which receives as input the measured value 
Y(k) and the input to the plant U (k) , and has as output the 
estimated value R(k) . This is multiplied by the measurement 
function H and the result is compared with the observed 
value Y(k). The difference is then multiplied by the gain 
G (k) and summed to the predicted value to produce the 
estimated value, as given by Equation (31) . 

When we have more than one state variable, an advantage 
of the Kalman filter not shown here is that we do not need 
to measure all state variables but only a smaller number 
to get all the estimated values we need to feed back to the 
plant . 

As an example, suppose we have values from Table I in 

/v 

that a = 0.7, b = 0.001, T = 100 and R q = 0.8; suppose also 

that the variance in the estimated value of R is M = 0.1, 

o 

and that the variances in the random input and in the measure- 
ment are known to be Q = 0.1 and S = 0.1. Let our measure- 
ment equation be given by 
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FIGURE 7. Flow graph in the time domain for the 
rework system with a Kalman Filter 



51 



Y(k) 



R(k) + v (k) 



which in terms of Equation (29) , implies H = 1. With these 
values we can compute the Kalman filter gain for four 
periods using Equations (32) to (35) : 

G (0) = 0.5, 

G (1) = 0.58, 

G (2) = 0.60, 

G (3) = 0.60, and 

G (4 ) = 0.60. 



Suppose the measured values are 

Y ( 1) = 0.75 , 

Y (2) = 0.82 , 

Y (3) = 0.86, and 

Y (4) = 0.85 . 

and we guess a value for the initial reliability 

A 

R(l|0) = 0.724. We obtain, by using Equations (30) and (31), 
the values 
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0.737 , 



R(l) 



R(2) 



0.825 



R(3) 



0.859 , and 



R(4) 



0.856 . 



Notice that the reliability before the rework, R(k) , 
is not defined for k = 0. For this reason we use as initial 



guess R(lj0) the value of the reliability before the first 
rework computed from the initial value R q = 0.8 and use the 

/N /N 

gain G(0) to estimate R(l), the gain G(l) to estimate R(2) 
and so on. The gain would not be delayed by one period if 
we had decided to call the reliability before the first 
rework by R(0) . 

D. EXTENSIONS TO MORE THAN ONE DIMENSION 

In a missile we can consider for example its booster, 
cruise motor, guidance, and warhead as four subsystems that 
are independently influenced by the environment and have 
different reliability deterioration factor b. To extend 
our model for this case we need to decompose our system 
in several distinct subsystems. Here we may put the 
equations for the subsystems in matrix form. 

We can define the vector of reliabilities of the 
subsystems that will be called our state variables , as 
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R(k) 



R ± (k) 
R 2 (k) 
R 3 (k) 
R 4 (k) 



and a vector of reference reliabilities as 



?D 



**02 

^3 

^D4 



The input vedtor can now be defined as before in terms 
of the effectiveness of the rework for each subsystem: 



i 

a 

£ 

i 




*D1 - R l' 


| u 2 (k) 


= [a^ a 2 a 3 a^] 


CM 

1 

CM 


1 

U 3 (k) 


^3 R 3 


U 4 (k) 




R D4 ” R 4 



U(k) = a T [R d - R(k) ] . 
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We can then write our state equation for the system with 
random input as 



/s 

R(k+1) = $ R (k) + <£ U(k) + w (k) , (36) 



where the matrix <j> is 



$ 




0 

0 

0 




0 



0 



0 




J 



The measurement equation is 



Y(k) = H R(k) + v(k) 



(37) 



and the vector of estimated values is given by 



R(k) = R (k | k) = R (k | k-1) +G(k)[Y(k) -HR(k|k-l)] 

(38) 



and the vector of predicted values by 



^ — "hT ^ -hT 

R(k | k-l) = e ~ R (k-1 | k-1) + e ~ U(k-l), 

(39) 



where G(k) is the vector of the Kalman filter gains. 
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If we are interested in the overall system reliability 
and if as in the missile case the subsystems are inde- 
pendent and in series, we have the output equation 

4 

R(k) = tt r (k) 
i=l 

which can be added to the control signal flow graph as 
shown in Figure 8. For subsystems arranged in parallel 
or series-parallel, the system can be handled in a similar 
manner, and needed changes in the output equation will not 
influence the system itself. 

In this chapter we applied tools from control theory 
to our reliability model and showed how to find the limits 
of the rework effectiveness a for which the system is 
stable and the final value of the item reliability when k 
becomes large. We also showed how a Kalman filter may be 
applied to smooth the estimated values of item reliability 
and feed back this estimated values upon which the decision 
about a is made. 

In the next chapter we will show how the decision 
about a can be optimized with respect to a given performance 
measure . 
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w(k) 



v(k) 




R(k) 



n 

& R-j_(k) 

i=l 



FIGURE 8. Signal flow graph for multidimensional rework 
control system 
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V. OPTIMIZATION 



In the previous chapter we have studied our control 
model for the case we have linear equations and a single 
level of rework effectiveness a. In this chapter we will 
extend this study for non-linear equations and several 
levels of rework effectiveness, and show how to find the 
sequence of rework that optimizes a given performance 
criteria. First we will discuss various performance 
criteria that could be used for our control model, and then 
we will show how to find the optimal sequence of rework 
levels by dynamic programming. 

A. THE PERFORMANCE MEASURE 

The optimal control problem in our thesis is to find 
a control u(k) as defined by the expression 

u ( k ) = a (k) (R D - R (k) ) , (40) 

which causes the system to follow a trajectory that opti- 
mizes a performance measure J. The rework effectiveness, 
a(k), is now allowed to change from step to step and in 
fact the problem now becomes one of finding the sequence 
of rework effectiveness a(k) that yields the optimal 
sequence of control u(k) . Notice that since a (k) is no 
longer constant. Expression (40) is non-linear. 
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First let us consider a trajectory (that is, the 
sequence of values the item reliability follows over time) 
that attains the desired reliability with a smallest num- 
ber of reworks. This type of problem is called a "Minimum- 
Time" problem and the performance measure to be minimized 
may be generally stated as 

fc f 

J = / dt = t f - t Q 

fc o 

for the case of our discrete system, this becomes 

N-l 

J = It, 
k=0 



where (N-l) is the first period the desired reliability, 

R_ , is attained. 

If instead we want to minimize the deviation of the 
first state of our system from its desired value we have 
the type of problem called a "Terminal Control" problem. 

Here, possible performance measures are 

N 

J - l R d - R si (k) ' 
i 

where the summation is done in all dimensions or reliabili- 
ties of the sub-systems. If positive and negative deviations 
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are equally undesirable these deviations should be 
squared . 



To transfer a system from an arbitrary initial state 
R q to some specified desired value with minimum expenditure 
of control effort or with minimum cost, we need to minimize, 



N-l 

J = l C (k) T 
k=0 



N-l 

l cu(k)T , 
k=0 



where we assume that the cost of a rework C(k) is propor- 
tional to the control u(k) or equivalently that the cost 
is proportional to the reliability improvement, R (k) - R(k), 
after each rework. 

Since these criteria are completely distinct in their 
concepts, we shall follow the usual approach of using a 
combination of them, that for the one dimensional problems 
takes the form, 

N-l 

J = HN[R d - R g (k)] + l [Q ' + cu (k) ] T 

k=0 

or 

N-l 

J = HN[R d - R g (k)] + l [Q + Gu (k) ] , (41) 

k=0 
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where HN, Q and G are relative weights so that by adjusting 
their values we can weight the relative importance of each 
criteria with respect to the other. These weighting fac- 
tors can be functions of time if the relative importance 
varies with time. Notice that now the performance measure 
is a combination of criteria and has no physical meaning. 

It should be noticed also that the rework effectiveness 
a(k) now has the definition given at the end of Chapter III 



R s (k) - R(k) 
a(k > ‘ ' R D - ROc) 



and that even when Rp is not equal to 1.0 we do not need 
the decision rule 



M 



0 if R D - R (k) < 0 

< 

a (k) if R - R(k) > 0 



because the constraints in a(k) do not permit a to become 
negative . 

B. OPTIMIZATION USING DYNAMIC PROGRAMMING 

We now wish to show how the performance of our rework 
control system may be optimized in terms of the combined 
performance measure (41) . An optimal solution would be 
expressed as a sequence of reworks that minimizes this 
criterion. We shall approach the problem with the method 
of dynamic programming, as developed by Bellman. 
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The basic notion here is given by Bellman's principle 
of optimality: 

"An optimal policy has the property that whatever the 
initial state and initial decision are, the remaining 
decisions must constitute an optimal policy with regard 
to the state resulting from the first decision" [5] . Thus 
if N (R(k+l)) is the minimum cost to reach the final 

state at k = N, starting from the state R(k+1) at time 
t = (k+l)T, and is the cost to move from kth to the 

(k+l)th state, then C£ N (R(k) ,u (k) ) , the minimum cost to 
go from the kth to the (k+l)th state when we use a particu- 
lar control u(k) is given by 



C* N (R(k) ,u(k) ) = J k/k+ 1 ( R (k) ,u(k) ) + J k+1/N (R(k + 1) ) • 

(42) 

The optimal decision at the instant k, u*(k) , is the 
decision that minimizes Cj* (R (k) , u (k) ) over the set of 
possible controls u(k). 



J* .(R(k) ) 
kN 



min fC* (R(k) ,u(k) ) ] 
u (k) 



(43) 



Expressions (42) and (43) form the functional equations 
of our dynamic programming approach. The optimal sequence 
of decisions are built up from the final state N backwards 
toward the earlier states. This is necessary in order 
that J* be known prior to the calculation of C*. The values 
of R(k) are given by our model, and the constraints of the 
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problem are 



0.0 < R(k) < 1.0 



and 



0.0 £ a(k) £ 1.0 

The first step in the computational procedure, then, 
is to find the optimal policy for the last stage of opera- 
tion. This is essentially a matter of trying all of the 
allowable control values at each of the allowable state 
values. To limit the required number of calculations, and 
make the computational procedure feasible, the allowable 
state and control values are discretized. The degree of 

I 

approximation depends on the separation of these discrete 
values and on the method of interpolation used and can, 
of course, be adjusted. 

For each discrete value of R(N-l) we try all discrete 
values of u(N-l) and calculate the resulting state R(N) . 

The optimal control for this rework is the one which yields 
the minimum cost. The procedure is repeated for all the 
other discrete values of R(N-l) . This gives a table of 
optimal policy for each value of R(N-l) at the last stage. 
Since the cost J , „ is dependent on the value of the 
state R(N-l) and on the value of the input applied, u(N-l)/ 
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the minimum cost J* and the optimal control u* (N-l) are 
dependent on the value of the state R(N-l) • For the last 
rework we use only the first term (terminal control) of 
our performance measure (41), but for the other, successive 
stages we must compute each term in that expression. When 
a state does not coincide with one of the discrete values 
we have to use interpolation to find the corresponding value 
of the performance measure. Since a direct search is used 
to solve the functional equations, the solution obtained 
is regarded as the global minimum. 

A flow chart describing the computational procedure is 
shown in Figure 9. A Fortran program for this flow chart 
is described in Appendix B. 

As an example of the use of this dynamic programming 
procedure, we sought the optimal sequence of reworks for 
the case where b = 0.001, T = 100, and constraints, 

0.8 £ R (k) <_ 1.0 

and 



0.7 < a (k) < 1.0. 



We chose to work with twenty discrete values for R(k) , 
thirty discrete values for a(k), and four reworks. We 
input in the program the equations of our model : 



j READ input data: nunber of 

' stages = N T ; nurrber of state 
! values - NOC; timber of 
control values = NTJ; other 
f ’-ecuired info oration 



L 



I K = 0 

Calculate and store j£-m 
i'utnirrun cost of last stage) 
for all discrete values of 
3'N) 



_L 



j K = K + 1 

Bet 300 equal to the lesser 
•discretized value by making 
•i = 1 



Set COSMTN to a large 
positive number; SET U(k) 
equal to the lesser 
discretized value by 



I 



! CAIjCJlATC the value of R(N-Ktl) . 

! Use this value to select the 
appropriate stored value of 
! ^ (interpolation required j 

‘ if hot' coincides with a discrete , 
| 'value) ; CO’ 2’jTC C* the vunirarn cost } 
j over the final Ol-K) stages fcr 
i 'J(M-K) 



Glance U(K)j to , 
j the cent 
discrete control 
. (J = ■>!) 




5T-RE GUSi-llN 
COSTOP (N’-K, 3) and 
CO; in UCPTOI-X,!) 




STORE the value C* in 
OOSiOOT ; STORE the value 
U* in EMDJ 



PRUTT optimal ( 

I controls j 

t EOPTC^-k.I) and 
‘ rdnimuni costs , 
j COSTOP (N-K, 3) for ! 

) all discrete points: 

| and all stages. 



a 



FIGURE 9. Flow chart of the computational procedure for 
the dynamic programming solution for optimal 
control of the rework control system 
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R (k) + a (k) (l-R(k) ) 



R s (k) = 

R (k+1) = R s (k) e " L , 

and considered as input variable the rework effectiveness 
a(k) since in this program the equations need not be linear 
and this is the variable of immediate interest for the user. 
For the performance measure J given by Expression (41) 
we considered four distinct types of problems obtained by 
changing the weights. These are: 

(1) Minimum Cost, so that the weights of the combined 

criterion J (41) are 

G = 1.0, Q = 0.0 and HN = 0 . 0 , 

(2) Terminal State, or minimizing deviations from the 
final value. Here the weights in J are 

G = 0.0, Q = 0 . 0 and HN = 1 . 0 , 

(3) Terminal State and Minimum Cost, where we assume 

a desired reliability R D = 0.99, and a relatively important 
terminal state represented by the weights 

G = 0.9, Q = 0 . 0 and HN = 500.0 , 
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and 



(4) The linear combination of the three criteria 
represented by the weights 

G = 0.05 , Q — 0.05 and HN = 500.0. 

The dynamic programming solutions for the four cases 
are presented for each rework period in the Computer Output 
Section of this thesis. From these solutions we can now 
find, for each of the four cases, the optimal sequence of 
reworks for a given initial reliability R q (entering with 
this value at the first rework table and, interpolating 
among the values listed, we can follow the process until 
the final state is reached) . 

The results are in Tables II to V for two values of 
initial reliability of an item. As we can see from these 
tables, the sequence of reworks that minimizes our perfor- 
mance measure depends on the value of the initial reliability 
R q , for given values of the desired reliability and of 
weightings in the performance measure. 

Summarizing, we have applied a dynamic programming 
procedure to our control model of the reliability of an 
item and were able, for the four selected cases, to find 
an optimal sequence of reworks that minimizes a selected 
performance measure. For other selections of performance 
measure J, relative weightings HN, Q and G, constraints, 
initial reliability R Q and model parameters b and T, we 
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TABLE II 



RELIABILITY FOLLOWING 


REWORK 


FOR INITIAL RELIABILITY 


V 


AND OPTIMAL SEQUENCE OF 


REWORK LEVELS a (k) 




Case (1) - 


Minimum 


Cost 








N-l 


PERFORMANCE MEASURE: 


J = [Rd 


-R (k) ]H + l Q+GU (k) 
k=0 


with 


H = 0.0 , Q = 


0.0 , 


G = 1.0 and U(k) = a (k) 




R =0.8 
o 






k 


R s (k) 


a (k) 


R s (k+1) 


0 


0.800 


0.60 


0.890 


1 


0.890 


0.60 


0.924 


2 


0.924 


0.60 


0.933 


3 


0.933 


0.60 


0.936 


4 


0.926 








R =0.9 
o 






k 


R s (k) 


a (k) 


R s (k+1) 


0 


0.900 


0.60 


0.926 


1 


0.926 


0.60 


0.933 


2 


0.933 


0.60 


0.936 


3 


0.936 


0.60 


0.937 


4 


0.937 
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TABLE III 

RELIABILITY FOLLOWING REWORK FOR INITIAL 
RELIABILITY R , AND OPTIMAL SEQUENCE OF 
REWORK LEVELS a (k) 

Case (2) - Terminal State 

N-l 

PERFORMANCE MEASURE: J = [R_-R_(k)]H + l Q+GU(k) 

13 k=0 

with R d =0.99, H = 1.0 , Q=0.0, G = 0.0 

and U ( k ) = a ( k ) 



R = 0.8 

o 



k 


R s (k) 


a (k) 


Rg (k+1) 


0 


0.800 


0.90 


0.972 


1 


0.972 


0.90 


0.988 


2 


0.988 


0.90 


0.990 


3 


0.990 


0.90 


0.990 


4 


0.990 









R =0.9 
o 




- 


k 


Rg(k) 


a (k) 


Rg(k+1) 


0 


0.900 


0.90 


0.982 


1 


0.982 


0.90 


0.989 


2 


0.989 


0.90 


0.990 


3 


0.990 


0.90 


0.990 


4 


0.990 







TABLE IV 



RELIABILITY FOLLOWING REWORK FOR INITIAL RELIABILITY 
R q AND OPTIMAL SEQUENCE OF REWORK LEVELS a(k) 

Case (3) - Terminal State with Minimum Cost 

N-l 

PERFORMANCE MEASURE: J = [R_-R (k)]H + l Q+GU(k) 

U k=0 

with Rj^ 0.99 , H = 500.0 , Q = 0.0 , G = 0.9 

and U(k) = a(k) 



R_ = 0.8 

o 



k 


R s (k > 


a (k) 


Rg(k+1) 


0 


O 

• 

00 


0.60 


0.890 


1 


0.890 


0.60 


0.924 


2 


0.924 


0.60 


0.933 


3 


0.933 


0.867 


0.979 


4 


0.979 









R 

o 


0.9 




k 


r s oo 


a (k) 


R s (k+1) 


0 


0.900 


0.60 


0.926 


1 


0.926 


0.60 


0.933 


2 


0.933 


0.60 


0.936 


3 


0.936 


0.864 


0 . 978 


4 


0.978 
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TABLE V 



RELIABILITY FOLLOWING REWORK FOR INITIAL RELIABILITY 
R , AND OPTIMAL SEQUENCE OF REWORK LEVELS a(k) 

Case (4) - Terminal State in Minimum Time with 
Minimal Cost 



N-l 

PERFORMANCE MEASURE: J = [R_-R (k)]H + £ Q+GU(k) , 

U S k=0 



with R d = 0.99 , H = 500.0 , Q = 0.05 , G = 0.05 

and U(k) = a(k) 



R = 0.8 

o 



k 


R s (k) 


a(k) 


R s (ktl 


0 


0.800 


0.60 


0.890 


1 


0.890 


0.60 


0.924 


2 


0.924 


0.868 


0.977 


3 


0.977 


0.90 


0.989 


4 


0.989 







R = 0.9 

o 



k 


R s (k) 


a(k) 


R s (k+1) 


0 


0.900 


0.60 


0.926 


1 


0.926 


0.60 


0.933 


2 


0.933 


0 . 835 


0.973 


3 


0.973 


0.90 


0 .988 


4 


0.988 
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would obtain different sequences of rework effectiveness 
a(k) that optimizes the selected performance measure. A 
suitable selection of the performance measure and relative 
weightings being a responsibility of the manager according 
to his purposes, experience and data available. Thus the 
cost of each rework, the minimum value allowable to the 
item reliability, minimum and maximum levels of rework that 
are physically feasible to implement, and other factors 
like the urgency to attain a desired reliability in minimum 
time, or a constraint in the budget, are taken into account 
when selecting the weightings in the combination of all 
the criteria. 

In the next chapter we will present the general conclu- 
sions of this thesis and suggest the areas where one could 
do further study involving the rework control model. 
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VI. CONCLUSIONS AND RECOMMENDATIONS FOR FURTHER STUDY 



A simple, closed- loop inventory storage system from 
which the stored items are removed, sent through a rework 
mechanism to improve their reliability, and returned to 
storage has been investigated in this thesis. Such a sys- 
tem might be a stock of large ordnance which is acquired, 
stored and reworked at regular intervals . 

A general rework model was developed by using the tools 
of control theory and a solution in closed form was found 
that permits further analysis of the system in a much easier 
way than if we had a solution in the form of series. This 
model also permitted us to find in a direct way the values 
of effectiveness of the rework that makes the system unstable 
and the final value of the item reliability over several 
reworks when the system is stable. 

A Kalman filter was used in the control model to obtain 
an estimation of the item reliability when we have random 
error in our measure of the item reliability or when the 
rework process introduces randomness in the item relia- 
bility. This model was extended to more than one dimension 
for the case where we have systems composed of several 
subsystems in series, parallel, or both. Finally a study 
was made for the case where we can have several levels of 
rework and want to find the optimal sequence of such reworks 
that minimizes a performance criteria that leads to the 
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case of minimum time, minimum cost, terminal state or a 
linear combination of these criteria. 

This thesis shows that control theory may be applied 
for administrative problems where a mathematical model des- 
cribing the system can be derived. Control theory gives 
a broader understanding of the problem and due to its flexi- 
bility would permit us to include in the model many other 
variables not studied in this thesis. 

Among the several interesting areas which might be 
pursued in further study, one is to determine methods for 
measuring or estimating the item reliability by frequent 
diagnostic checks so that this measurement process can be 
incorporated in the model. Another area is to broaden the 
field of the rework model so that other state and input 
variables could be incorporated. These include the proba- 
bility that an item chosen at random might have reliability 
exceeding some predetermined reliability requirement, and 
other input variables (such as items expended for training 
purposes, use, or obsolescence, as well as new items acquired 
periodicaly) . The rate of rework, the rate of expenditure, 
or the rate of acquisition could then be sought as control 
variables besides the rework effectiveness. Non-linear 
functions caused by the decision rule of doing a rework 
only when the item reliability is smaller than the desired 
value, can also be solved in terms of non-linear control 
theory . 
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A reliability control model has been developed in this 
thesis for a system of stored items requiring rework. 

It is hoped that the results presented here will not only 
be useful to inventory managers and high-level planners 
but will also generate further interest in the application 
of control theory to administrative problems. 
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APPENDIX A 



Z TRANSFORM 

A. 1 DEFINITIONS AND PROPERTIES 

The solution of difference equations by the z-transform 
method is very useful, because we can transform difference 
equations into algebraic equations in z. Once solved 
this algebraic equation we can then find the inverse z 
transform of this solution to obtain the solution in the 
time domain. 

Thus, given a discrete function of t, f(kT), the 
z-transform of this function, symbolized by F(z), is defined 
by 



Z [f (kT) ] = F(z) = l f (kT) z k (44) 

k=0 



From this definition we can obtain several useful 
properties. For example, let's find the z-transform of 
f ( (k+1) T) : 



Z [f (k+l)T) ] 



l 



f ( (k+1) T) z 



k=0 



£ f ( kT ) z k+ 



k=l 



= z[ l f (kT) z~ k - f (0) ] 
k=0 



z F(z) - zf ( 0 ) . 
-(45) 
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Thus when a difference equation is transformed into an 
algebraic equation in z by the z-transform method, the 
initial data are automatically included in the algebraic 
representation . 

Other useful properties are 

z [af (kT) ] = a F (z) , A 



and 



z [l(kT) ] - , ( 47 ) 

where l(t) is the unit-step function: 

for t < 0 

for t > 0 

These two properties permit us to find the z transform 
of a constant, since any constant a can be expressed by 
the product al(kT): 

z [a] = az/z-1 (48 ) 

v Another property is the distribution property 
z[f 1 (kT) + f 2 (kT) ] = z[f 1 (kT) ]+z[f 2 (kT) ] (49) 



0 

1 (kT) = < 

1 
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The z-transform of the most common functions are given 
in several references listed in the Bibliography [1] , [4] 

and [10 ] . 



The inverse transform of F(z) is f (kT) , and is denoted 

by 



From the tables of z-transform we can find the inverse 
of the simplest function. More complicated z-transf orms 
may have to be expanded into partial fractions so that the 
tables can be used. Normally, we expand F(z)/z instead of 
F ( z ) into partial fractions because this leads to functions 
with z in the numerator after we multiply back by z, and the 
functions of z appearing in tables of z transforms usually 
have the factor z in their numerators. 

When the partial fraction expansion does not give 
tabulated functions, we may have to find the appropriate 
z-transform from the definition. This will be the case 
in our problem. 

Similar properties of z-transforms exist for inverse 
z transforms and will be useful: 



Z 1 [F (z ) ] 



f (kT) 



(50) 



Z _1 [F(z)/z] 



Z 1 [z ^F (z) ] 



f [ (k-l)T] 



(51) 



Z _1 [aF(z) ] 



aZ -1 [F(z) ] 



af (kT) 



(52) 




a l(kT) 



a 



for t > 0 , 



(53) 



and 



Z _1 [F 1 (z) + F 2 (z) ] 



Z _1 [F 1 (z) + Z _1 [F 2 (z) ] . 



(54) 



The characteristic equation of a system is equal to 
the denominator of the overall transfer function and the 
stability of the system can be determined from the location 
of the roots of the characteristic equation [10] . A con- 
dition for stability is that all roots must lie inside the 
unit circle or 



Another important property of the z transform is the 
Final Value Theorem [10] : 

"If x(t) has the z transform X(z) and X(z) has no 
poles (roots) outside the unit circle (|z| <1, that is the 
condition for stability) , then the final value of x(t) 




(55) 



or x(k) is given by 



lim x(t) 



lim x(k) 

k-«» 



lim[ (z-l)X(z) ] . " (56) 

z-*-l 



To prove this note that 



oo 



Z [x(k) ] = X (z) = l x(k) z 

k=0 



79 



Z[x(k+1) ] = ZX(z) 



Zx ( 0 ) 



l x(k+l) 
k=0 



z 



-k 



Hence 

zX ( z ) - zx(0) - X ( z ) = (z-l)X(z) - zx(0) 

CO 00 

= j. x(k+l) z ^ £ x(k) z ^ 

k=0 k=0 

from which we obtain 



(z-l)X(z) = z x(0 ) + l [x(k+l) - x(k) 3 z 

k=0 



Because of the assumed stability condition, we obtain, 
as z -*■ 1, 



lim [(z-l)X(z)] = x( 0 ) + x(oo) - x(0) = x(«) 

z->-l 

which is Equation (56). 

A. 2. Z TRANSFORMS OF EQUATIONS (8) AND (9) 

We found that the reliability immediately before the 
(k+1) th rework is given by equation ( 57) . 

-bT -bT 

R (k+1) = (l-a)e R(k) + ae (57) 
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— bT -bT 

Since (l-a)e and aR^e are constants, we can apply 
the properties (45), (46), (47), (48) and (49) to find the 

z-transform of both sides: 

— bT — bT 

Z [R ( (k+1) T) ] = Z [ (1-a) e R(kT)] + Z [ae ] 

— bT — bT 

= Z [R ( (k+1) T) ] = (1-a) e i Z[R(kT)] + ae Di Z[l(kT)] 

zR(z) - ZR(1) = (l-a)e bT R(z)+ ae bT , (58) 

z~l 

which is Equation (12) of Chapter I when we use R(l) as 
the initial condition. 

Similarly for Equation (8) : 

-bT 

R s ( (k+1) T) = (l-a)R s (kT)e + a (59) 

zR (z) - zR(0) = (l-a)e bT R (z) + a y (60) 

s S Z _L 

which is Equation (10) of Chapter I. 

A. 3. EXPANSION INTO PARTIAL FRACTIONS 

Equations (58) and (60) can be solved for R s (z) 
and R(z), giving: 



R s (z) = 



z- (1-a) e 



-bT 



■R(0) 



az 

[z- (1-a) e bT ] [z-1] 



( 61 ) 
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and 



-bT 

ae z 



R(z) -bT R ^ + -bT 

z- (1-a) e DX [z- (1-a) e Di ][z-1] 



(62) 



In order to have simpler expressions we can proceed now to 
expand into partial fractions. We will exemplify it for the 
last equation only, since they differ only by a constant 
factor in the last term. 

The first term has already a simple form so that we 
need only to expand the last term. For reasons stated in 
Section A. 2, we will expand this term divided by z: 



-bT 

ae z 



[z- (1-a) e ] [z-1] 



ae 



-bT 



[z- (1-a) e bT ] [z-1] 



A 

z-1 



+ 



B 

z- (1-a) 



-bT 



e 



(63) 



To determine the coefficient A we can multiply both 
sides by (z-1) and evaluate the expression at z = 1 : 



ae 



-bT 



z- (1-a) e 



-bT 



= A + [- 



B ( z-1) 



z=l 



z- ( 1-a ) e 



-bT 



z=l 



This gives 



A 



1 - 



-bT 



ae 

(1-a) 



e 



-bT ‘ 
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— bT 

To evaluate B we multiply both sides by [z-(l-a)e ] 

-bT 

and evaluate at z = (l-a)e , so that the term with A 
becomes zero: 



ae 



-bT 



z-1 



z= (1-a) e 



-bT 



0 + B , 



this gives 



B 



ae 

(1-a) 




-1 



/ 



Substituting back into (63) , we have 



ae 



-bT 



ae 



-bT 



(z-l) [z- (1-a) e -bT ] [l-(l-a)e- bT ][z-l] 



ae 



-bT 



-bT -bT ' 

[ ( 1-a) e -1] [z- (1-a) e ] 



and R(z) can now be put in the form: 



R ( z ) 



z- (1-a) e 



-bT 



•R (1) + 



ae 



-bT 



z-l ,, ,, , -bT, 

(l-(l-a)e ) 



-bT 

z ae 

[z-(l-a)e J [(l-a)e -1] 



(64) 



The expansion for R s (z) follows immediately from this last 
equation: 
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z 




z- (1-a) e 



z 



a 



( 65 ) 



[z- (1-a) e 




A. 4. INVERSE Z TRANSFORM 

The first and last terms of Equations (64) and (65) 
are not found in the tables of z-transform so that we will 
proceed now to develop the expressions of these z-transforms 
from a guessed form of f (kT) . 

Since the first and last terms are identical in terms 
in z, differing only by constant factors, we need to do the 
development just for the expression in z. 

Thus, let's assume that 



f (kT) 




Applying the definition of z transform we have 



F(z) 



l f(kT)z _k = l (l-a) k e 
k=0 k=0 



k -bkT -k 
e z 



( 66 ) 



Multiplying both sides by (l-a)e 



-bT -1 



, we have : 



z 



, , > -bT -1„ . . 

( 1-a ) e z F (z ) 
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Subtracting (67) from (66) side by side we get 



[l-(l-a)e bT z _1 ]F(z) 



l (l-a) k e" bk V k - l (l-a) k+1 e- b(k+1)T z 
k=0 k=0 



or 



[l-(l-a)e~ bT z 1 ]F(z) =1+1 (l-a) k+1 e“ b(k+1)T z~ (1+1) 

k=0 

- I (l-a) k+1 e- b(k+1)T z" (k+1) 
k=0 



or 



[1- (l-a) e bT z -1 ]F(z) = 1 



Thus 



F(z) = 



, , -bT -1 

1- ( l-a) e z 



z- ( l-a) e 



-bT 



( 68 ) 



From this it follows that 



,-l 



z- (l-a) e 



-bT 



= f (kT) = 



,, . k -bkT 

(l-a) e 



(69) 



Also , 



z- 1 t 



a 



. -bT ... . -bT . 

z-(l-a)e (l-a)e -1 



] = 



,, . k -bkT 

a(l-a) e 

, , , -bT n 

(l-a)e -1 



(70) 



- (k+1) 
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and we can get similar expressions by changing the constant 
part. From (52), (54), (69) and ' (70) we can now find 

the inverse z-transform of equations (64) and (65): 



R(kT) = (1-a) k e“ bkT R(l) + 



ae 



-bT 



1- (1-a) e 



-bT 



,k -b (k+1) T 
a(l-a) e 

, . . “bT m 

(l-a)e -1 



or 



R(kT) = (l-a) k e" bkT R(l) + — 



-bT 



-bT 

(1-a) e -1 



r , . . k -bkT . , 

((1-a) e -1] 



(71) 



And from (65) 



R (kT) = (1-a) k e bkT R + 
s o 



a 



-bT 

(1-a) e -1 



r , k -bkT , . 

[(1-a) e -1] (72) 



These are the equations (13) and (11) of Chapter I, respec- 
tively, and constitute the solution of the difference 
equations (57) and (59) in the time domain. 
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APPENDIX B 



COMPUTER PROGRAM DISCUSSION 



In Chapter V we have shown that the functional equa- 
tions of the dynamic programming to compute the optimal 
costs over the set of possible controls are 



The program that performs these computations is 
listed in the Computer Program Section of this thesis and 
is composed of three subroutines and two functions . The 
subroutine QUANTU computes the discrete values of the state 
variable R(k) , and of the inputs U (k) . For each discrete 
value, the subroutine STATE computes the value of the state 
variable at the next state and the dynamic equations of 
our model (15) and (16) has to be input in this subroutine. 
With the vlaue of R(k+1) , the subroutine cost selects the 
value of the performance measure J by interpolation between 
the two nearest discrete values stored in the array of 
discrete values of optimal cost, COSTOP, of the (k+1) th 
state. This is possible because the computation is done 



C* N (R(k) ,u(k) ) 



J k,ktl (R(k) ' u(k) ) +J k+l,N (R(k+1)) 



and 




min [C* .(R(k) ,u(k) ) ] 
u (k) XN 
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backwards from the last stage of operation. The costs of 
the last stage are computed by the Function HN, and the 
subsequent cost for each state and each value of the input 
U are computed by the Function G. These costs are summed by 
the subroutine COST and in the main program the minimum 
value is found for all discrete values of U(k). These com- 
putations are repeated for all the discrete values of the 
state variable as can be seen in the flow chart of figure 9. 

The program output presents the optimal costs and corres- 
ponding effectiveness of rework a to go to the final state 
for all discrete values of the reliability so that starting 
with a given initial reliability at the initial state or 
first rework table we can, interpolating through the values 
of the next state until the final state, find out the optimal 
sequence of rework for that particular initial condition. 

The program output for the examples of Chapter V are 
presented in the Computer Output Section of this thesis. 
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COMPUTER OUTPUT 



CASE 1 



9 A T A IN 


PUT 




N NX 


f\" ) 


MODc 


+ 20 


10 


J 


N XM AX 


UN II 


•JMA/ 


3 1 . 0 


o.s 





costs )F m’unfi a ' 1 o ] p y i a 



FABLE r JR " r M E =l\AL STATE 



FINAL 


v I N I '• J 1 


STATE 


COST 


i • o'j 


0 .0 


■J. 9 9 


0. 0 


J . c 6 


J . J 


0.97 


0 .0 


0. 9b 


0. 3 


0 . ) 'j 


^ . 0 


0. 9 A 


0.0 


0. c 3 


\ n 
0.0 


J . -7 


c.o 


0 . 9 » 


0. 0 


J. 39 


c.o 


0 . 33 


0 .0 


0. 57 


J. 0 


0. 3o 


0.0 


0 . -33 


0.0 


3 . 3-A 


0.0 


0. a 3 


0. 0 


0. 32 


0.0 


0 . 3 l 


0.0 


0.30 


0. 0 
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JT'CJL LAW 





COMPOTE 


9 OUT 0 JT 


CASE l 




7 A 3 L E FOB TnE 


4 T H 0 EU 3 '< X 




CURRENT 
R S ( X ) 


N : XT 
R3 ( Xt i) 


MI Mi MUM 
C 33 T 


OPT i ML 
RE WORK 


1.00 


0. 962 


0.4 


0 . 60 


C. 93 


0. 9 32 


G.-f 


J.SJ 


0.9 9 


J • 9 54 


0.4 


C. 60 


1 » 
1 

r- i 
c* 1 

• | 
o t 
1 
1 
1 


0 . ? 5 1 


0. + 


G. 60 


0.9 c 


0.947 


0.4- 


0.60 


0.93 


0 . 9 92 


0.^ 


o ♦ C J 


vj • 9 -+ 


G • 9 39 


0. •+ 


C . 60 


0.99 


C. r 35 


■J .-t 


0.6 0 


G • 3 2 


0.921 


0 • 4 


v • J J 


0.9 1 


0 • ^ j. 3 


3 • •■+ 


G • £ J 


yJ • 'J ~ 


U .9 34 


J • "-h 


G • 5 J 


0.3 3 


V-> • ■* /„ vj 


0 • -t 


C . 6 J 


0. 3 7 


J.-ij 


O.-r 


0 ♦ SO 


0 . 3 o 


0.9 12 


0.4 


C .cO 


G . o 


0.909 


0.- 


0 . 60 


0 • 3 v 


0# 3 G5 


0.3 


G • 5 0 


G. Si 


0 .9 01 


0.9 


0 .60 


0 . 3 ’ 


0 • 3 3 7 


0 .+ 


0 • D 0 


,-\ -« 1 
O • O _1» 


0 . 3 9 3 


0 . + 


0.60 


o q "i 

^ • -> 'J 


0. 390 


0 . 9 


w.oO 
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C^PIJTER OUTPUT 



CASE 1 





T A 3 Lt FOR THE 


2TH RERQRis 




clr.re rr 

R3 ( K ) 


NeXT 
*■% o ( K +■ ± ) 


M MI V IM 

COST 


optical 

RE WGRrt 


1. CC 


0.9 62 


C . 7 


G. 60 


0.05 


0.5 53 


0. 7 


0 . 60 


C. 9 0 


0.904 


C . 7 


C.5 5 


0.97 


0.951 


0. 7 


C . 6 0 


C • *56 


0. ? 97 


0.7 


G. SO 


0.9 5 


0. o 43 


0. 7 


0.60 


C • 9 -i 


0. 5 25 


C. ^ 


G. 3 J 


o : 


0 . 5 3 5 


0.7 


C • 60 


C . 9 1 


0. 5 2 _ 


C. 7 


0.60 


0.9 1 


0.52 3 


0. 7 


G . 60 


L • £ 'z 


0 • ~ 3 '+ 


0.7 


C. .60 


0 • So 


C. ^30 


0. 7 


0 . 60 


c. c 7 


G • 0 l o 


0.7 


C . 6 0 


0. 06 


0 . *> 1 2 


C. 7 


U . 60 


C. 85 


U * 7 0"* 


0.7 


0 . 60 


G* 3 4 


0. 9C5 


0. 7 


0 . 60 


G » o 0 


0. 5 01 


0. 7 


0 . 60 


G. 3.2 


0 • 3° 7 


G. 7 


0 • 60 


C. 81 


0 . 3 53 


0. 7 


0.60 


0. 3 0 


0. R 9 0 


0. 7 


0 . 60 
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CC 4, ?'JT ER 'j'JTDijr 



CASE 1 



TAiLE FC’ ’HE 2 T H RE-HRX 



C URR E 1 T 
R S ( K ) 


NEXT 
■ 7 S U+l) 


'•‘I NI NUN 

COST 


op r i ^ al 
RE NO 2k, 


* .0 0 


0.982 


I. i 


0. 2 0 


C. 99 


0.2 52 


1.1 


C.60 


2 


J . - 5 f 


l. i 


C. 80 


0. c 7 


0.-51 


1.1 


C. 20 


0.96 


0.947 


i. • — 


C. 80 


G. 9 5 


w . ^ 93 


1 . i 


C . Oo 


0. - A 


0 • 7^7 


1. A 


0. 6v> 


0.9 2 


0.935 


l.i 


C . 50 


0 . 9 2 


0 .9 ii 


1. 1 


0. cO 


0.9 i 


0. -28 


1 . 1 


C. 50 


0 • ? 


0 • ? 2 V 


1. A 


0.60 


0 . 6 5 


0 . ? 20 


l.A 


C . dO 


0."7 


0 . c 18 


i • 4 


C . C \J 


0. 36 


0.912 


l.A 


O * 3 J 


0 • fi 3 


0 .9 09 


1. 1 


0 ♦ CU 


0.S4 


0.905 


i.i 


C • 0 0 


o . a 2 


0.901 


l.A 


< •> 
Vj • C vJ 


0 . 82 


0. 597 


l.i 


C. 50 


3.31 


0.392 


1 « X 


C. 50 


0.80 


0. 5 90 


1. 1 


G • 6 0 
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COMPUTER T JT D U T 



CASE i 





TA3 15 =CR ’H = 


IT H 33..<0?K 




CURRENT 


M EX T 


«rn 'i 


OPT IV ?L 


p 5 ( K ) 


X S( K+l> 


COST 


9 E W C P is 


1.03 


0. 3c? 


1. t 


0. j0 


C. 9 9 


0. 3 55 


i . 


0.60 


0. cs 


0.953 


l. + 


0.50 


C. 9 7 


0 . 9 5 \ 


• , 


0.60 


r> r-. ' 

O • - 


u • 9 4 7 


1.5 


0 . 60 


C • 9 3 


0 .9 A 5 


1. ? 


0.60 


,•> r> *■ 

U# 7 -*■ 


0.93? 


1.4 


0 • 60 


0.9 3 


J . n 3 5 


l. V 


0.63 


o " 
^ • 7 


0. 931 


1.5 


0. 50 


C. 91 


0. 9 35 


1.5 


0 .59 


0« C J 


0. ^ 24 


1. + 


0 . 60 


0. 3 3 


0 • -* 1 j 


t • T 


0.60 


0. 6 T 


0.3 15 


l. 3 


0. 50 


C. 86 


0.913 


1. + 


0 . 50 


G ♦ 8 3 


0 .909 


i • 4 


0.60 


0.3 4 


0.Q05 


i.-* 


0 .60 


3.9) 


0 • c 0 1 


i 


0.60 


3. 3 3 


0.337 


1 . 4 


0 . 60 


C. ? 1 


0. ^ -33 


1*4 


0.60 


C . 3 3 


0.390 


1 .4 


0.60 
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C TAUTER OUTPUT CASE 2 



DATA INPUT 



N 


NX 


NU 


MODE 


4 


lO 


30 


0 


X AIN 


XM4X !J M I N 


U M 'AX 


U • S 


i.o 


0.6 


0.9 



COSTS IF DERATION A ID OPTMAl CC'JTR 1L LAW 
TAdlE FOP ThE FINAL S^ATE 



FI'IM J IN I MUM 

STATE COST 



i .00 


3 . 0 J 


0.9 9 


0 .00 


0 . ?3 


0.00 


0 . 


0 .00 


0 • 96 


0 .0 0 


0 . 95 


0.00 


0.^4 


0 . 0 u 


0 . 93 


0 . 00 


0.92 


L . 01 


0.91 


0.0 1 


0.69 


0 . 0 '. 


0.63 


0.01 


0 . 3 1 


0 . 01 


0 . 66 


0.02 


0.35 


0.02 


0 • 3 ■+ 


0 .0 2 


“)• 33 


0.03 


0 . 33 . 


0.05 


J . ^ i 


0.0 3 


J . 30 


0 . Oh - 
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COMPUTER OUTPUT 



CASE 2 



TA5LE FOR THE 4 T H REWORK 



CUR* z IT 


NEXT 


MI MUH 


OPT l > AL 




P 5( K ) 


RS (K+l) 


COST 


RE U0RK 




1.00 


0 . 9 39 


0 . 0 


0.59 




C. 99 


0.9 90 


0. u 


o • ^ 0 




0.93 


0.959 


0. J 


0 . 90 




C. 97 


9. 9 33 


0. 9 


0.30 




0. 9c 


0. 9 S 7 


0. J 


0.90 




C. 9 5 


•J. 9 36 


0.0 


0 .90 




0. 9 j 


0.^53 


0. 0 


3 . 30 




C . 9 3 


0. 9 34 


*j . 0 


0 . 90 




Q. 92 


0. 3 33 


0. 0 


0.2Q 




0.91 


0. 9 32 


0. 0 


0 . 9 } 




1 

1 

1 

1 o 
1 • 

( sJJ 

1 L 
1 
1 
i 


0.331 


■J • J 


0.90 




U • 3 3 


C. 9 30 


0.0 


0.3 0 




0.3 7 


0.3 79 


0. J 


0.2 0 




0. 5 6 


0 . 273 


0. 0 


0.9 0 




0.3 5 


0. 5 77 


9. J 


0.20 




0. 3 4 


C • 9 76 


0.0 


0 • 3 J 




1 

1 

1 

1 o 
1 • 

1 iX 
1 l • 

1 

1 

1 


0.975 


0 • J 


u • w J 




0. 97 


0 .9 74 


0.0 


0.20 




0. 31 


0 • 9 73 


0 .0 


0.90 




■> O > 

o» Uu 


0.9 7 2 


0.0 


0.90 
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COMPUTER OUTPUT CASE 2 



TA3LE FOP THE Z T H R E//7R K 



C CRR ? IT 
R S ( K ) 


NEXT 

S ( K * i ) 


NT NI MU 7 
CAST 


JPT i v AL 
7 E WORK 


i. • 0 0 


C • 9 90 


0. J 


0.90 


C. 99 


o • 9 9 0 


0.0 


0. =0 


0.93 


0.959 


0.0 


0.90 


C. 9 7 


0. 3 58 


0.0 


C.90 


0. 9 6 


0. 9 c7 


0. J 


G . 30 


C. 95 


0 • ^86 


0.0 


0 • 7 J 


0. 9 7 


0 • 3 c 5 


0. J 


0 . -0 


C . 9 3 


0. 9 34 


0. 0 


G . 90 


0. 9 * 


0 • 9 c 3 


0* 0 


0 . 90 


G. 9 J. 


0. r >32 


0. j 


G.9J 


C. 3^ 


0.9 57 


0. J 


0 . 90 


0. S3 


C. 9-0 


0. 0 


0 • 90 


0. 3? 


0.979 


0.0 


0.P0 


0. 3 6 


0 • 9^3 


0*0 


G . 3Q 


0. 3 5 


0*977 


0. 0 


u . 90 


C • 5 ** 


0. 976 


0. 0 


0 . 30 


0* 83 


0.9^5 


0. 0 


u . 90 


u • d c. 


0 • 9 7 * 


0. J 


0.30 


0. Si 


0.9 73 


0. J 


0* 90 


0. 3 C 


0 . 9 72 


0. 0 


•u . -0 
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COMPJTER OUTPUT 



CASE 2 





TA3LE FCR THE 


27 H REWORK 




CURRENT 
RSI M 


NEX T 
RS (K, +i ) 


M I M I mu 3 
COST 


OPTIMAL 

R t W OR K 


1 . DC 


o.?-o 


0.0 


0.90 


0.9 9 


0.9 90 


0.0 


G. 30 


0.9 3 


0.9 39 


0.0 


C. 5*J 


0.9 7 


0.9 33 


0 • vj 


o*90 


0 • 9 6 


0.937 


0.0 


0.90 


0.9 3 


0 • j 86 


0. 0 


C. 90 


J . 99 


0.9 35 


0.0 


0.90 


0.93 


0.-39 


0. 0 


C. 50 


i 

1 

O' 1 

• 1 
O 1 
1 
1 
1 


0.9-33 


0.0 


C. 90 


C. 51 


G . 9 32 


c. 0 


0.50 


0 • 8 ^ 


0.9 31 


0. 0 


0 . 5 J 


C. 9 • 


0 . 9 S 0 


p ) 

VJ • 


0.90 


0.5 7 


0.979 


0.0 


C . 50 


C. Si 


0.973 


0. 0 


C.90 


0.35 


0. 977 


0.0 


C. 90 


C. c'-+ 


0. P?3 


u * 0 


C.90 


0.33 


0.975 


0.0 


C.90 


0 • C 3 


0. 974 


0 .0 


0. 90 


0 . 8 1 


0.9 79 


0. 0 


0 . 30 


w • E c 


0. 972 


0.0 


G. 90 
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CC‘-1P'JT£P TUI D Ur 



CASE 



T A 3 L t F CP THE 1TH RE'OU 



CUR?, E IT 


NEXT 


M 1 1 1 1 3'J3 




OPTIMAL 


RS(K ) 


9 S< K+i) 


c OS 7 




REWORK 


1 

1 

1 

1 

| • 

1 o 

1 C< 

1 

1 

1 


0.990 


0.0 




0.90 


0.9 9 


0.990 


0 .0 




0.30 


C. 9 i 


0 • 3 S 3 


0.0 




0 .30 


C. 9 T 


0 • 9 33 


9 .0 




0.90 


vj • ' sc 


0.937 


0.0 




0.^2 


0.9 3 


0.9 36 


0 .0 




C. 50 


C • 5 + 


0.2 55 


0.0 




0 .50 


0.93 


0. ?34 


0 . 3 




C .30 


C.92 


0.933 


0 . 0 




0 .°Q 


3 . 9 i 


J • i 3 2 


0 .0 




G .90 


G • S 7 


0.-31 


0.0 




C • 90 


0.8 3 


C . 7 30 


0 .0 




0 . 7 J 


0.3 7 


0.579 


0.0 




Ls . 7 J 


0.36 


0.973 


0.0 




C. 90 


0.65 


0.977 


0.0 




0.90 


0.34 


0 • 9 76 


0. 0 




C . 7 0 


0.3 3 


0. 9 75 


0 . 0 




0.90 


0 . 8 2 


0.974 


0. 0 




G • 9 0 


G • c i 


0 . '73 


0 . J 




C. 90 


o .ro 


0.9 72 


0.0 




C. 9 0 
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COMPUTER OUTPUT 



^ ^ 3 £ 





OAT A If 


IPI.JT 






N 


MX 


MU 


MICE 




4 


2 J 


CO 


0 


A. MI 


N 


XMAX 


UM I i'l 


UM aX 


0 « 


a 


1.0 


0.6 


0.0 



J F )PE r A~ 


I )M A ID OP T i '1 A 4_ 


1 -* 6 L z c C R 


~H2 F I i\AL 3 T A T : 


F INAL 


‘U MI ■’ JM 


ST A fE 


C 35 T 


1.00 


0.05 


0.-9 


J .00 


0. 93 


O.Oo 


o. or 


J. 25 


0. 96 


0.52 


0. 95 


o.<u 


O.o'+ 


j. . -vi 


0.92 


1 . 0 5 


0.9 2 


2.75 


0 . ' l 


3.5-3 


0.3 2 


V • 3 4 


J. 33 


5.50 


J.. 37 


s5 . 7 o 


0.35 


^.04 


0. 35 


c • 4*+ 


0.9 V 


10.9 4 


0 . 5 3 


i - • ^ 5 


0 . 52 


19. 27 


J . 3 A 


16.11 


0.30 


1 3.05 



C N T R J L L 4 W 
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i ki 



C IMPUTE 0 OUTPUT 



i c 



5 





T43LE F0° TmE 


4TH R EW 13 X 






CURREUT 
R 5 ( K ) 


NEXT 
R S ( K.+ 1 ) 


V I Ml MUM 
C 1ST 


OP 

Rc 


T 1 ML 
'aORX 


1 .00 


u . 3 7 7 


0 . D 




0 .73 


C. 9° 


0.9 79 


C . o 




0.30 


0. 9 J 


0 .9 7 3 


0.7 




0.32 


c. = / 


0.-/9 


0.7 




j . 3 o 


1 ' 
O . ^ O 


0 • ^ / 5 


0.7 




0.3- 


0.95 


C. - 73 


0.7 




<j . i ^ 


G . r + 


o . - r = 


0. / 




0.35 


0 .p 3 


0.--- 7 9 


0 . 7 




^ -* v 

U • w J 


C . 9 0 


0.^/9 


0. 0 




C . 33 


0 . c 1 


0. ■» 30 


0 • d 






0.39 


0.77" 


0.3 




0 . j - 


0 .3 3 - 


0 • 7 £ 0 


0 . 3 




C. 30 


0.87 


0. p ’9 


0.3 




0 . 90 


0,3 3 


0 . P 73 


0. 3 




^ . "5 '«y 


0. £5 


0.9 77 


0.8 




'w . ? 'sj 


0.84 


0. P75 


0.3 




C • 9 0 


0.53 


0.975 


0.3 




C • 30 


0 .3 2 


0.9 74 


C. 9 




C . ^ G 


0.51 


0. 9 73 


0.9 




0. 90 


0.30 


0.972 


C. 9 




0. 90 



100 



COMPUTER 1J TP JT CASE 3 
TABLE FOR THE 3 T H REWORK 



CURRE NT 
K 3 ( l\ ) 


HEX T 
3 S( K+l) 


Mi -i r-vjM 
CIST 


J 0 T I v A i_ 
5 E a ] rl K 




1. OC 


C. 9 6? 


1.0 


0 • D J 




C. 9 9 


0.953 


1.0 


0 . 5 7 




C.9 3 


0.35 + 


l.u 


0 , 50 




0. 9 7 


0.9 31. 


l. o 


j • 60 




0.96 


0.3+7 


U 0 


J .63 




L • ^ o 


0. 9 42 


l . 0 


G . 9 3 




G • 9 9 


0. 3 3° 


1. 1 


J • 5 3 




> r- ", 

O. 7 -* 


0.9 35 


i • L 


j • ~ 




C. 9 2 


0. 3 3i 


1. 1 


0.r>3 




vJ. 9 I 


0. 9 33 


1 . 1. 


J . 9 3 




o ^ 

u • 0 ? 


0. 9 £3 


i • 


J .30 




0.33 


G • 9 70 


1. i 


0.5 3 




C. 3 / 


0.916 


l .1 


0.33 




o 3 i. 

Vj • o O 


o.3 1: 


i • i 


O • O 'J 




0. 9 5 


0.9 09 


1 .1 


3 • 60 




0. 3 t 


0.^05 


1 .1 


C.SJ 




O'? 
J • O 


0.501 


1 .1 


0.50 




0.3 3 


0. 3 3 r 


1 .1 


0.60 




G • 3 L 


G • 5 9 ~ 


T 1 

J • J 


0 .60 




3. 30 


0. 3 30 


i .1 


0.30 
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COMPUTES OUTPUT 



CAS c Z 



TA3LE F UP THE 2TH 9E.MR.< 



CURRENT 


■NEXT 


M I N I MUM 


OPTICAL 


r< s ( k ) 


RS ( K + 1 ) 


COST 


RE a/ JR X 


l.OJ 


0.PE2 


1.4 


J. 60 


0.99 


J • 5 5 3 


1 • 


\J • - J 


C. 9 c 


0.9 34 


1.4 


J . 60 


0.97 


0.951 


1 . 4 


o • 3 \j 


C. 5 6 


C. °4T 


1 • -+ 


u « o ) 


C. 9 3 


0 • 9 A 3 


i • + 


0 . t 9 


Cm 9 9 


0 . 39 


i. V 


0.60 


0. 9? 


0 • 9 3 5 


1. 4 


J • 59 


G. 92 


0.901 


1. + 


J • 5 J 


C. 9 1 


O'. 5 38 


1. t 


0 . 0 J 


c* s 9 


0. 9 24 


l* 9 


0.80 


0* 8 3 


Cm 930 


i • T 


0.60 


0. S 7 


0 • 9 i e 


i. + 


U t w J 


G. SS 


0.512 


1 . v 


0.60 


C • 8 5 


0. 9 0? 


t. 4 


0 • 6 J 


C. 8 + 


0.905 


1. 4 


0. 50 


C. 82 


0. 901 


1.- 


0 • 6 J 


0. 8 ’ 


0 . 8 97 


1. 4 


J . 50 


C. 81 


0 . 3 93 


l.H- 


0 • 50 


C • 3 0 


0.3 90 


1. 4 


0.5) 
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COMPUTER OUTPUT CASE 3 



T46 LE FOP THE 17 H REWORK 



C UP-1 E 1 T 


NEXT 


MI hi mum 


OPT I v AL 


RS (K ) 


RS U+i) 


COST 


R E R 0 R K 


i .GO 


0 . - 33 


i. r 


C. 5 0 


0.99 


0 . 9 58 


l. r 


C . 60 


0.9 3 


0 . -5- 


1, 7 


C. 6 0 


• 7 t 


0.931 


1. 7 


0.6 J 


0 . c 5 


0."-,' 


1. 7 


o • C J 


C • 9 5 


G • ■ + j 


1.7 


C.-'j 


0.9 + 


0.939 


1. 7 


C. 6 0 


i 

1 

i 

1 o 

i * 

1 vO 
1 u> 

1 

i 

1 


•> ~~cr 
O • 7 2 J 


1. 7 


L. • D <J 


0.? 3 


0. c 21 


1. 7 


C • c 0 


0.91 


0 . 4 C 3 


i. 7 


0 .60 


1 C ' 
J • *-j V 


G . 9 


I. 7 


0. 60 


0.3c 


C.9 20 


1.7 


o • 6 J 


J.S7 


0.9 16 


:. 7 


G. 30 


0*56 


0.911 


1. 7 


C . zO 


C. 33 


0.909 


1 . 7 


2 • 3 j 


0.34 


J . 9 G 5 


1 • 7 


G .60 


C. S3 


G • "T 1 Li J> 


1 . 7 


G • 6 0 


0.32 


0. 397 


t 7 

v • i 


G . 60 


G • £ i 


0 • 3 9 3 


1. 7 


0.6 0 


J . 3 G 


m. 390 


L. 7 


C . cO 
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COMPUTER OUTPUT ' ;s = a 
OAT4 INPUT 





NX 


NU 


MODE 


A 


20 


30 


U 


X M l M 


X M AX 


U Vl I N 


'l'- 4 AX 


0.3 


1.0 


C.6 


0.9 



0 . 



0. 
0 . 
0 . 



0 . 



3 "RATION 


A NO 'jPTIMAL 


FOR T h E 


final state 


•A L 
i T5 


MINIMUM 

COS’ 


► 00 


0.05 


► 9 9 


0.00 


, 9 9 


0.06 


► ^ 7 


0.23 


• _ 9 


0.52 


*7 2 


0. oi 


9 4 


1 . 4_ 


93 






i • Or* 


L. 


1 . 7 5 




5.59 


.■> *? 


4.3- r 


c ° 
^ o 


-.60 


37 


6. 7o 


36 


S .04 


35 


= .4t 


34 


1C. 94 


3 5 


1 ?.. 55 


£ '? 


14.07 




16. 1 i 


30 


1^.05 
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C T'l PIJ T E R TJTPJT 



CASE 4 



TA8LE F OK THE 4 T H PEOOXK 



CU~2E 


NEXT 


*1 MI V(JM 


OPTI M,4L 


KS( k ) 


XS(X^i) 


CO S T 


9E0G2K 


1. GO 


0.9 29 


0.1 


1 

0.39 


0.39 


0.9 90 


0 . 1 


C . 9 0 


C. 9° 


0.939 


n •) 

vj * L 


C. 90 


0.97 


0. 9^3 


o. i 


3.90 


C « > d 


0.957 


o.i 


7.30 


0.9 5 


0.9E6 


0. i 


0 ^ ( J 


C • 9 3 


0 . 9 35 


0.1 


. 17 -J 


J • 9 3 


U . i 


0.1 


C . 90 


C.92 


0.983 


0 . 1 


C. 90 


0 • 9 i 


0 . 9 22 


0. 1 


0.90 


0. 3 9 


0.9 31 


3.L 


0 . 90 


J • 3 3 


C. 9 80 


0.1 


0.90 


0.8 7 


0. 9 79 


0. 2 


C. 90 


0.8 6 


0 . 9 7 3 


0.2 


C . 9 3 


0.3 5 


0.97 7 


0. 2 


<u • 3 0 


0. 34 


0 . 9 76 


0.2 


C. 90 


0.82 


0.375 


0 . 2 


C. 90 


0.32 


0.974 


0.2 


0.93 


0. 31 


0. 973 


0.2 


0.90 


0. 30 


0.372 


0.3 


0 .90 
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COMPUTER OUTPUT CASE 4 
TABLE FOR THE 3 r H REWORK 



CURRENT 
KS( K ) 


NEXT 
RS (K + l ) 


MINI MUM 
COST 


U P T I ^ AL 
REWORK 




l.OC 


0.962 


0.2 


u • oO 




0.9 9 


0.95 3 


0. 2 


C • 50 




0.93 


U.95 + 


0.2 


G • 60 




0.9/ 


0. 557 


0 .2 


0 • 7) 5 




C • 9 d 


0 .9 5- 


0.2 


C • 69 




‘J • 7* J2 


0. Ro- 


0 ._ 


C.7-5 




C • 9 •+ 


0.569 


0.2 


C .30 




0.5 5 


0. j^o 


0.2 


C • c5 




C* 9 0 


0.3 79 


0.2 


0 . 33 




J. Q L 


U • O 0 


0 .2 


0. 39 




vj • O * 


o.R n 


0.3 


0.33 




J . 2 1 


C. 5- 0 


) .2 


u • 7 0 




0.8 / 


J.3?9 


3.2 


0.90 




0. 3 5 


0.9 73 


0.2 


0.90 




0. 33 


0 . 9 77 


0.2 


0 . ^3 




0.8 4 


0. 


0 .2 


0.9 0 




0. 35 


0 • r 75 


0.2 


0.90 




0. 3 l 


0. 375 


0 .3 


0 . 90 




C. 81 


0.5 77 


0 .2 


0 .90 




0.2 0 


0.973 


0 .2 


0. 90 
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COMPUTER 0 UT° JT 



C a 5 E 4 



T A3 LE FOR THE 2TH REWORK 



CURREN r 


NEXT 


VI II VUH 


optical 


P S ( K J 


P. 5 ( K 1 ) 


COST 


R E W 0 k X 


1. OC 


0. 9 60 


0.2 


0.60 


C. 99 


0.9 53 


0.2 


0*50 


C. 93 


3*954 


C.2 


0.50 


C . 9 7 


0.951 


0. 2 


0 • c 0 


1 

1 

O I 

U' 1 

• 1 

r -J 1 

1 

1 

1 


0. 947 


0.2 


0 • 4 0 


0.9 5 


0. 943 


0. 3 


C • t> 0 


C. 94 


0.9 39 


0.3 


0 • CD 0 


n - 

U • 7 _ 


0.955 


0. 3 


G. 60 


C. 9" 


0 • 331 


0.3 


0*50 


0.9 1 


3*923 


0. 3 


G . 2 : J 


C *3 ^ 


0*324 


0.3 


0.60 


0. 3 3 


0.9 20 


0.3 


0 . 60 


0. 37 


0.916 


0. 3 


0*50 


0# 3 6 


0. 911 


O -I 
'J • - 


0 .50 


0.35 


0. 3 09 


0.3 


0.60 


0. 3 4 


0*905 


0.3 


0.60 


0# 3 3 


0. 3 01 


0. 3 


0*60 


0. S i 


0.557 


G. 3 


u * j0 


0.81 


0. 3 93 


0. 3 


0 • o 0 


1 

1 

-Tj | 

• 1 

o \ 


0. 3 90 


0.3 


0 * so 
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C CP PUT 


ER OUTPJT C 


ASE 4 




7A3LE F CP, TH 


E 1TH REWORK 




CURRENT 
P. 3 ( M 


H EX T 
RS(K-t-l) 


MINI M ' ) M 
COST 


OPT I MmL 
RE WORK 


1.00 


0. 9 62 


0.3 


G . 60 


C • 9 9 


0. 958 


u. 3 


0.60 


C. 93 


0.9 59 


0. 5 


0.60 


C. = 7 


0.951 


0. 2 


. 0 0 


0.9 c 


0.9^7 


0. 3 


D. 20 


C. 93 


0. =42 


0. 3 


0 • 6 J 


C. 9- 


0 . 9 29 


0. 3 


C. 60 


0.9 3 


0. 975 


0. 3 


0.60 


G . 9 2 


0.931 


0.3 


G . 60 


C. 91 


0.923 


0.5 


0 • 50 


0 .39 


0.924 


0.3 


0 . ^0 


C. d 3 


0 . - Z J 


0 . 7 . 


0.60 


0. 3 7 


u • 913 


0. 5 


0 . c J 


0.38 


0.9 12 


C. 3 


0.60 


0.3 5 


0.909 


0.3 


C. 60 


0.3*+ 


0.905 


C. 3 


J . oO 


0.33 


0.9 01 


0. 5 


u . i 3 


C. 3 2 


0. 397 


0.3 


L) . 0 0 


0.81 


0.3 95 


0.3 


C • 60 


C . 3 C 


0.390 


G. 3 


0.63 
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COMPUTER PROGRAM 
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OOOOOOOOOOOOOOOOOUOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 
O' rn >t*in o r- x ch o *-4<\i ^ or- ao Ch o ^ cm <r» in o r- ao o o 1-9 <\j m ^ in or- oo & o»-«cm m *4" in x 
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oooooooooooooooooooooooooooooooooooooooooooooooo 
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